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ABSTRACT 

We present XMM-Newton observations of a complete sample of five archetypal young radio- 
loud AGN, also known as Compact Symmetric Objects (CSO) or Gigahertz Peaked Spectrum 
(GPS) sources. They are among the brightest and best studied GPS/CSO sources in the sky, 
with radio powers in the range L5GHz=10 43 ~ 44 erg s _1 and with four sources having mea- 
sured kinematic ages of 570 to 3000 years. All five sources are detected, and have 2-10 keV 
luminosities ranging from 0.5 to 4.8x 10 44 erg s _1 . A detailed analysis was performed, com- 
paring the X-ray luminosities and Nji absorption column densities of the GPS/CSO galaxies 
with their optical and radio properties, and with those of the general population of radio galax- 
ies. We find that, 

1) GPS/CSO galaxies show a wide range in absorption column densities with a distribution not 
different from that of the general population of radio galaxies. We therefore find no evidence 
that GPS/CSO galaxies could reside in a significantly more dense circumnuclear environment 
such that they could be "frustrated" radio sources — hampered in their development. 

2) The ratio of radio to X-ray luminosity is significantly higher for GPS/CSO sources than for 
classical radio sources. This is consistent with an evolution scenario in which young, compact 
radio sources are more efficient radio emitters than large extended objects, at a constant ac- 
cretion power. 

3) Taking the X-ray luminosity of radio sources as a measure of their ionisation power, we 
find that GPS/CSO sources are significantly underluminous in their [OIII]^^ line luminos- 
ity, including a weak trend with age. This is consistent with the fact that the Stromgren sphere 
should still be expanding in these young objects. If true, this would mean that here we are 
witnessing the birth of the narrow line region of radio-loud AGN. 

Key words: galaxies: active - X-ray: galaxies 



1 INTRODUCTION 

Ever since their discovery, it has been speculated that those com- 
pact radio sources that show convex-shaped radio spectra at cm 
wavelengths, may be young objects fchklovskvl 1 96llBlakel 1 970). 
As a class, they were named Gigahertz Peaked Spectrum (GPS) 
sources after their characteristic radio spectrum (see O'Dea 1998, 
for a review), most likely caused by synchrotron self absorption 
iFanti et aljfl990l ISnellen et alj EoOO, e.g.). High resolution Very 
Long Baseline Interferometry (VLBI) observations have shown 
that these sources are typically up to a few hundred parsec in 
size, often exhibiting jet and/or lobe structures on two opposite 
sides from their central core — the reason wh y they are also 
called Compact Symmetric Objects (CSO; eg. Wilkinson et al. 
Il994l) . The most compelling evidence that GPS/CSO are indeed 

* E-mail: j.vink@astro.uu.nl 



young radio sources comes from VLBI monitoring observations, 
showing that the bright archetypal objects i n this class have hot- 
spot propagation velocities of ^0.1 -0.2c lOwsianik & Conwavl 
Il998tl0wsianik et alJl99llTschager et all2000l) . indicating kine- 
matic ages of ~10 3 years. This is in contrast to speculations that 
GPS/CSO galaxies are small due to confinement by a particularly 
dense and clumpy interstellar medium (ISM) that impedes the out- 
ward propagation of the jets Ivan Breugel et alll984tlO'Dea et all 
Il99lt) . 

Note that a large fraction of GPS sources, in particular those 
showing a convex radio spectrum peaking at higher frequencies 
than a few GHz, turn out to be identified with high redshift quasars 
(z~2-3). Their connection with the population of GPS/CSO galax- 
ies at lower redshift is not clear, and they may well be a completely 
separate class of object which just also happen to exhibit a convex 
shaped spectrum ISnellen et alll999l) . By no means it has been es- 
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Table 1. The sample of GPS/CSO radio sources with |6| > 20° from the Pearson & Readhead 1 1988) catalogue. Indicated are, (column 1) the coordinates, 
(column 2) redshift, z, (column 3) radio flux density at 5 GHz, Srq hz (erg s —1 ), (column 4) radio luminosity at 5 GHz, Lsghz (erg s — 1 ), (column 
5) the 5007 A[OIII] emission line luminosity, £[oiii] ( er g s_1 ) from lLawrence et alJ 1 19961). and (column 6) the kinematic age of the radio hot spots 
IPolatidis & Conwavl200i except for B1358+624) 



Source 


Position 




2 




log ^5GHz 


lo g£[OIII] 


Kinematic Age 




(J2000) 






Jy 






yr 


B0108+388 


0l' l ll m 37.3 s 


+39°06'28" 


0.668 


1.6 


44.0 


40.8 


570 ± 50 


B0710+439 


07 h 13 m 38.1 s 


+43° 49' 17" 


0.518 


1.6 


43.7 


42.4 


930 ± 100 


B1031+567 


10 h 35 m 07.0 s 


+56°28'47" 


0.45 


1.3 


43.4 


41.6 


1800 ± 600 


B1358+624 


14 /i 00 m 28 gs 


+62°10'39" 


0.431 


1.8 


43.5 


41.8 


2400 ± 1000 b 


B2352+495 


23 h 55 m 09.4 s 


+49°50'08" 


0.238 


1.5 


43.0 


41.3 


3000 ± 750 



a Obtained from the NASA/IPAC Extragalactic Database (NED) 

b The kinematic a ge of B 1358+624 has not been directly measured, but is based on on its size and the average size-age relation of GPS/CSO sources in 
IPolatidis & Conwa\ji2003T) . 



tablished that the GPS quasars may also represent a young stage of 
radio source evolution. 

Here we present XMM-Newton X-ray observations of a 
small but complete sample of all five GPS sources from the 
IPearson & Readhead ll988l) sample with Galactic latitude \b\ > 
20° (TableG). These are among the brightest GPS/CSO galaxies in 
the sky. All these sources appear to be young radio loud AGN with 
four out of five sources having measured kinematic ages of their 
hot spots, indicating ages of up to ~ 3000 yr. No kinematic age 
measurement for B1358+624 exists, but base d on its size and the 
age-s ize measurements of GPS/CSO sources I Polatid is'& Conwavl 
2003) we estimate its age to be 2400 yr. Although the sample is 
limited in size, it allows us to study the correlations between their 
X-ray, optical and radio properties, and to assess how they com- 
pare to those of mature radio galaxies. One expects that the X-ray 
luminosity is largely a manifestation of the instantaneous accre- 
tion power of the central black hole, whereas the radio luminosity 
is expected to evolve s ubstantial ly over the life time of the source 
(Read head et alJll996tlSnellen et alJl200d e.g.). Moreover, X-ray 
absorption measurements are an excellent means to probe the cir- 
cum nuclear density of the GPS galaxies. 

Several papers on X-ray observations of GPS sources have ap- 
peared in the literature, but one should be cautious to interpret these 
results in terms of X-ray pro perties of young ra dio-loud AGN. A 
first success was obtained by lO'Dea et all2 000) with ASCA. Al- 
though they did not detect B2352+495, they obtained a firm de- 
tection of GPS galaxy B1345+125 (PKS1345+125). Furthermore, 
iGuainazzi etal] 120041) detected B 1404+288 (Mkn 668). 1 Both are 
low redshift GPS galaxies exhibiting strong optical line emission 
and are powerful infrared emitters, and may be not representa- 
tive to the class of young radio-loud AGN (yet no reliable ages 
have been measured for these sources). The combined X-ray and 
radio observatio ns o f the G PS quasar B0738+313 presented by 
ISiemiginowska et alJ 1200 3l) show that, with its prominent kpc- 
scale X-ray/radio jet, it is certainly not a young radio-loud AGN. 

In order to easily compare our results for GPS/CSO galax- 
ies with the X - ray p roperties of a large sample of AGN by 
Sambr una et alJ 1 19991) we adapt here a cosmology with Ho = 
75 km s -1 Mpc -1 and q = 0.5. 

1 Prior to the submission of this publication we learned that Guainazzi et 
al. have detected a number of other GPS/CSO galaxies in X-rays, which do 
not overlap with our sample i Guai nazzi et all2005t) 



Table 2. Log of the XMM-Newton observations. The MOS exposure time 
and event rates refer to the average value for MOS1 and MOS2. 



Source 


Observation ID 


Start Date 


Exposures 


Event rates 








(MOS/PN) 


(MOS/PN) 






d/m/y 


ks 


ct s -1 


B0108+388 


0202520101 


09/01/2004 


16.4/12.0 


1.6/13.6 


B07 10+439 


0202520201 


22/01/2004 


14.2/11.2 


3.2/26.6 


B1031+567 


0202520301 


21/10/2004 


22.2/12.8 


34.4/93.0 


B 135 8+624 


0202520401 


14/04/2004 


12.4/12.0 


23.2/120.6 


B2352+495 


0202520501 


25/12/2004 


15.8/12.8 


3.7/28.5 



NOTE - The event rates refer to the total detector count rates, not the 



source count rates. 

2 OBSERVATIONS, SPECTRAL ANALYSIS AND 
RESULTS 

XMM-Newton Hansen et alJ 1200 ll) observed the five GPS/CSO 
sources as part of its guest observation program from January to 
December 2004 (Table [2j. All observations were made with the 
"Thinl" optical blocking filter. For the data reduction we used 
the standard XMM-Newton software package SAS v6.0.0. Unfor- 
tunately several observations were plagued by a high particle back- 
ground (see Table |2}. In the case of B0108+388, B0710+439, 
B2325+495 we removed time intervals with a high background 
count rate, using cut off rates of 15.5 cts^ 1 and 2.5 cts -1 for 
resp. PN and MOS. For B 103 1+567 and B 1358+624 the high back- 
ground persisted throughout the observation, and we simply used 
all available data. The observation B2325+495 had an intermediate 
background activity, so we selected time intervals with < 40 ct s _1 
and < 5 cts -1 for PN and MOS. For spectral extraction we used 
circular extraction regions with radii of 15" for B 103 1+567 and 
B1358+624, and 25" for the other three sources. Background spec- 
tra were obtained from rectangular regions near the source position, 
but excluding regions around 35" of the source. We extracted spec- 
tra for the two MOS CCD cam eras jTurneret aTll200il) . and the 
PN camera lStriide retall200ll) . For each source we combined the 
spectra of MOS 1 and MOS2, into one spectrum, which we anal- 
ysed using averaged instrumental response matrices. The potential 
systematic error introduced is small compared to the statistical er- 
rors, given the fact that the MOS1 and MOS2 are virtually identical 
instruments with similar instrumental response functions. 

All the five sources of the sample are detected. For the spec- 



The X-ray properties of young radio -loud AGN 3 



> X 

8 m 



™ o 

(0 T— 



CM 

X o 

CM 
I 



? B01 08+388 



H 1 1 — I — |- 



_1_ 



0.5 1 
B1 031 +567 



i o 



° T 
" o 



X o 

CM 



_1_ 



I 

0.5 1 
B2352+495 



-t-t-rf 



0.5 



1 



-i 1 r- 



_ B07 10+439 
c> 



X o 




H 1 1 — I — \ 



Energy (keV) 



10 0.5 1 

B1 358+624 



10 



Energy (keV) 



> °- 
% o 



CI) in 
c 1 

^ o 



o 
o 



H 1 — I — I — |- 



r 1 — tr- 1 - 



10 0.5 



10 



Energy (keV) 



Energy (keV) 



10 



Energy (keV) 

Figure 1. Observed XMM-Newton count rate spectra. The MOS spectra are shown in red and the PN spectra in black. 



tral analysis we employed a simple model consisting of a power 
law continuum and two absorption components: One represents the 
Galactic absorpti on, with an absorption column, JVh, fixed to the 
Galactic value of iDickev & Lockmanl Jl990h 2 . The other absorp- 
tion component corresponds to the absorption column intrinsic to 
the host galaxy. The redshift of this component was fixed to that of 
the galaxy, but the absorption column density was a free parame- 



2 We extracted the absorption columns from the online "./Vpi-tool" 

|http : / /heasarc ■ gsf c ■ nasa ■ gov/Tools| 



ter. The spectral analysis was done with the spectral fitting program 
xspec jAmauoj|l99r}) . using the absorption models of I Wilms et alJ 
( 2000) (called tbabs and ztbabs in xspec). 



The slope of the power law continuum was a free parameter 
for the high signal to noise spectra of B07 10+439 and B 1358+624, 
but fixed to 1.75 for the other three sources whose spectra are sta- 
tistically more limited. The best fit parameters, together with the 
inferred intrinsic X-ray luminosities between 2-10 keV are listed in 
Table|3| The spectra and best fit models are shown in Fig.Q 
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Table 3. Observational properties and parameters obtained by modeling the observed X-ray spectra in the range 0.5-10 keV. 





B0108+388 


B0710+439 


B1031+567 


B1358+624 


B2352+495 


PN 0.5-10 keV source count rate (10 3 ctss 1 ) 


4.7 ±0.9 


89.5 ± 3.0 


12.6 ±2.0 


44.3 ±2.7 


8.6 ± 1.2 


MOS1+2 0.5-10 keV source count rate (10~ 3 cts s -1 ) 


0.71 ±0.25 


33.0 ± 1.1 


4.3 ±0.5 


15.8 ±0.9 


3.4 ±0.4 


Galactic N H (10 20 cm- 2 ) 


5.80 


8.11 


0.56 


1.96 


12.4 


Normalization a (10 _5 ph s _1 keV _1 cm _2 @ 1 keV) 


3.3 ± l.l 6 


8.1 ±0.6 


1.2 ±0.2 


6.1 ± 1.6 


1.1 ±0.2 


Power law slope c (r) 


(1.75) 


1.59 ± 0.06 


(1.75) 


1.24 ±0.17 


(1.75) 


Intrinsic Nh (10 22 cm -2 ) 


57 ± 20 d 


0.44 ±0.08 


0.50 ±0.18 


3.0 ±0.7 


0.66 ±0.27 


Flux (2-10 keV) e (10^ 13 erg s^cm" 2 ) 


0.50 


4.0 


0.51 


4.8 


0.41 


Luminosity (2-10 keV)/ (10 44 erg s" 1 ) 


1.18 


2.16 


0.22 


1.67 


0.046 


C-statistic/bins 


137.5/99 


497.4/401 


104.2/97 


122.9/99 


126.3/99 



a Statistical errors correspond to AC = 1 (68% confidence limits). 

^ The 3a lower limit is 1.1 X 10 — 5 ph s _1 keV — 1 cm -2 . The source is detected at the 7a level. 
c Brackets indicate that the power law slope was fixed to this value. 
d The 3a lower limit is 1.8 X 10 23 era" 2 . 
e Including absorption. 

f Calculated for rest frame energies, ignoring absorption. 



44 - 



CO 

£? 

01 

X 42 



3 



40 - 



t — i — i — i — I — i — i — i — i — — i — i — i — i — — i — i — i — i — I — i — i — i — i — I — i — r 



T 

• BLRGs 
O NLRGs 
•B- RGs 
■ GPSs 



o o 



o ; 



J_£J I I I I I I I I I I I I I I I I I I I I I I I I L 



40 41 42 43 44 

Log Ls-iokev (erg s" 1 ) 



45 



44 - 



» 42 
S 

01 



eo 

3 



40 - 



38 



~ — I — I — I — | — I — I — I — I — | — I — I — I — I — | — I — I — I — I — | — I — I — I — I — | — r 

O /' - 



• BLRGs 

O NLRGs 

■0- RGs 

■ GPSs 



O 
* O • 

.6 • 



°o/t 



J I I I I I I I I I I I I I I I I I I I I I I I I I L 



40 41 42 43 44 

Log Lj.,0,,., (erg s" 1 ) 



45 



Figure 2. Left: Correlation between X-ray luminosity and radio luminosity of GPS/CSO sources (filled squares), compared to the radio-loud sample of active 
galactic nuclei bv lSambr una^^] ^9991). Righ t: Idem, for the [O III] line emission. In both figures the dotted lines are the average correlations for the radio 
loud sample obtained by Sambruna et al. 1 1999). The galaxies from the current sample are shown in blue, whereas Mkn 668 and PKS1345+125 are shown in 
green and red, respectively. . 



3 INTERPRETATION 

As we are interested in whether GPS/CSO galaxies are different 
from other radio-loud AGN we compare their X-ray properties 
to those of the sample of radio-loud AGN who se X-ray proper- 
ties were determined by ISambruna et all ^999) from ASCA ob- 
servations. We use all the sources of their broad line radio galax- 
ies (BLRG), narror line radio gal axies (NLRG) and radio galaxies 
(RG) subsamples. Since lSambruna et alJl 19991) do not quote upper 



limits for the non-detected intrinsic X-ray asborption components, 
we estimate upper limits proportional to the observed flux taking 
into account the errors on the detected instrinsic absorption com- 
ponents. 

In Fig. [2] we have set out the radio and [OIII] luminosities 
of the galaxies in our sam ple to their X-ray luminosities, and we 
compare them to the Sambruna et al. 1 1999) sample. In Fig.^we 
show the intrinsic absorption column density distribution. 
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Log N H 

Figure 3. The distribution of the intrinsic X-ray column densities for var- 
ious radio loud AGN. The lower panel are the r esults of the present study, 
the other panels have been taken from Sambruna et al. 1 1999). Upper limits 
are indicated by arrows. 



Table 4. A comparison of the X-ray derived absorpti on column J Vh and 
radio absorption column measurements A^hi IPihlstrom et all20O3ft . 





logAf H 


log A/hi 


B0108+388 


23.8 


21.9 


B 103 1+567 


21.7 


< 20.1 


B1358+624 


22.5 


20.3 


B2352+495 


21.8 


20.5 



These figures reveal three important properties of our sample 
of GPS/CSO galaxies: 1) for their X-ray emission GPS/CSO galax- 
ies are relatively radio-loud; 2) their [O III] emission is relatively 
low; 3) the column density distribution i s simi lar to those of radio- 
loud AGN classified by Sambrun a et aljfl99St) as narrow line radio 
galaxies (NLRGs) and radio galaxies (RGs), but the absorption is 
on average higher than those of broad line radio galaxies. 

As we discuss below these three properties support the hy- 
pothesis that GPS/CSO galaxies are indeed young radio-galaxies. 
Note that there are strong indications that GPS radio galaxies have 
relatively low [OIII] lumininosities with respect to their radio lu- 
minositi es as compared to compact steep spectrum (CSS) sources 
lQ'Dedll998h . This is probably the same trend between radio, 
[OIII], and X-ray luminosity that we report here, but without the 
X-ray luminosity as intermediary quantity. As noted in the intro- 
duction, the X-ray luminosity is the quantity that is probably the 
best indicator for the intrinsic power of the AGN. 

3.1 The absorbing column density 

An alternative explanation for the small extent of the radio jets in 
GPS/CSO galaxies that is still often considered in the literature 
JO'Deall998l) is that the radio jets are quenched by a high density in 
the vicinity of the nucleus, in other words they are "frustrated radio 
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Figure 4. The neutral hydrogen absorption column density, A^hi , toward 
the radio source IPihlstrom et alj2003t) versus the intrinsic X-ray absorption 
column density, A^h, for the four sources for which A^hi is measured. The 
solid line indicates A^h/x— ray = ^HI> me arrow indicates an upper limit. 



sources". It is clear from Tableland Fig.|3|that this is unlikely to be 
the case, as the intrinsic X-ray absorption is similar to other radio- 
loud AGN, with column densities ranging from a relatively modest 
4 x 10 cm to a considerably large 6 x 10 cm . A similar 
conclusion was drawn by Pihlstr 5m et alJ ^20 03") based on HI radio 
absorption observations of a large sample of compact radio sources, 
and by l O'Dea et al. 2005) from upper limits to the molecular gas 
content in GPS sources. Note that the X-ray absorption gives more 
stringent constraints on the actual gas column densities than HI and 
molecular absorption densities, as X-ray absorption depends on the 
total column toward the central source, whereas radio observations 
only probe the neutral fraction of the gas. This is quite an important 
distinction since AGN are expected to create an extended ionised 
region, as we discuss in section l3~2l Furthermore, X-ray absorp- 
tion probes the gas toward the accretion disk, whereas the radio 
absorption probes the neutral gas toward the radio source, which 
is situated at larger radii. It is therefore not surprising that for the 
four galaxies in our sample for which also HI absorption measure- 
ments have been made the HI column is always one to two orders 
of magnitude lower than the X-ray absorption column (Tableland 
Fig-HJ. Attributing the difference solely to ionisation effects would 
mean ionisation fractions of 90% to 99%. However, the ionisation 
fractions are likely to be lower, because a substantial part of the 
X-ray absorption may occur inside the central 100 pc, which is not 
probed by absorption toward the radio hot spots. 

IPihlstrom et aljj20 03) found a strong anti-correlation between 
the linear size of the radio emission and the HI column density, 
which they use to probe the average density profile of the inter- 
stellar medium. They do not consider ionisation effects, whereas 
this could be an additional cause for the observed anti-correlation: 
small jets are associated with young AGN, which do therefore not 
yet have an extended narrow line emission region of ionised gas 
(see section l3~2l . If this is the case it makes it less straightforward 
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to derive an average interstellar medium density profile from the 
relation between Nm and linear size of the radio emission, since 
the growth of the emission line region also depends on the density 
and UV luminosity of the AGN. 



3.2 The optical line emission: the case for an expanding 
emission line region 

Usually a high [O III] luminosity is taken as an indication for the 
presence of a powerful AGN, but Fig. indicates that GPS/CSO 
galaxies are relatively underluminous in [O III] compared to their 
radio or X-ray flux. This may not be too surprising if one takes 
into account that these are young radio galaxies, in which the AGN 
has switched on only a few thousand years ago. The reason is that 
it takes time to establish a large emission line region by photo- 
ionisation. This is best illustrated by a simple calculation, for which 
we assume an average interstellar medium density of 1 cm -3 and 
a typical luminosities of log L^q jjjj = 42, and log Lx = 44. The 

number rate of ionising photons (i.e. > 13.6 eV) is Nuv ~ 10 54 
ph s _1 for a power law spectrum with photon index -1.75. Com- 
paring this to the total nu mber of hydrogen atoms within a typi- 
cal region of 5 kpc radius ( Baum &Heckmarl 19891) . one finds that 
~ 10 67 atoms have to be ionised. In other words the central source 
has to shine for at least ~ 10 67 /10 54 = 10 13 s, or ~300,000 yr be- 
fore it has completely ionized a region with a radius of 5 kpc. This 
means that if an AGN has become only recently active it must be 
surrounded by a small, but rapidly expanding ionisation nebula (see 
IWhite et alJ2 003. for a calculation concerning the first generation 
of quasars). Hence, this would imply that in GPS/CSO galaxies we 
see the birth of the narrow line region of radio-loud AGN. 

In order to see whether this is the reason that GPS/CSO galax- 
ies are relatively underluminous in [O III] we consider a simplified 
model for the evolution of a Stromgen sphere. For an old ionisation 
nebula in equilibrium, the number rate of ionising photons (Nuv) 
should equal the number of recombinations, from which follows 
the equilibrium radius, Ri, of a Stromgen sphere: 



Ri — 



3AW 



(1) 



with an the hydrogen recombination coefficient. The initial rapid 
expan sion of the ionisation nebula is described by e.g. ISpitzeJ 



rf = Rf{l - exp(-n e a H t)}, 



(2) 



with n the radius of the ionisation front and t the time since the 
central source switched on. For small n e ant we can approximate 
this by 



r'l £s n a antR 3 



3N UV 



47rriH 

The [O III] line emissivity is then given by 



(3) 



N, 



[OIII] 
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(4) 



with / • the probability of emitting a photon at 5007 A after 
recombination. 

Assuming the interstellar medium densities in the GPS/CSO 
galaxies are more or less similar, we can expect the following cor- 
relation between the [O III] luminosity and the number rate of ion- 
ising photons for young GPS/CSO sources of kinematic age r, pro- 



vided that the age of the radio galaxy coincides with the birth of the 
ionisation nebula: 



-[OIII] 



oc tNuv ■ 



(5) 



As a first approximation we consider whether there is a rela- 



tion between the observables L 



[OIII] 



and r and Lx- However, 



Fig.|5|illustrates that the five sources in our sample do not support 
a simple scaling of £[QTJI] oc tLx- There may be various reasons 
why this is not the case: e.g. the interstellar medium density varies 
from galaxy to galaxy, their is no simple proportionality between 
Lx and Nuv, due to different spectral energy distributions (dif- 
ferent spectral slopes, or spectral breaks), and shocks induced by 
jet cloud interactions may provide an additional source of ionisa- 
tion. 3 Moreover, one of the outliers is B 1358+624, for which we 
only have an approximate age. 

However, a hint of what may the prime reason for devia- 
tions from Eq.[3]is provided by the very low [OIII] luminosity of 
B0108+388, because this is also the source with the highest ab- 
sorption column (Table|3}- So the most likely reason that the opti- 
cal emission is lower than expected is that the ionising UV flux is 
blocked by absorbing material close to the nucleus. The absorbing 
material is probably not the result of neutral hydrogen and helium, 
as, being close to the nucleus it would be ionised almost imme- 
diately, but dust grains, which may survive the extreme conditions 
close to the nucleus for 1000 yr to 10 6 yr, depen ding on the destruc- 
tion m echanisms and dust particle sizes (e.g. IVillar-Martm et all 
2001). For a young source like B0108+388 this means that an ap- 
preciable amount of dust may still enshroud the nucleus, frustrating 
the formation of a narrow emission line region. Note that this may 
also explain th e relatively high neut ral hydrogen column density 
of B0108+388 iPihlstrom et all2003l Fig.|4j; the hydrogen ionisa- 
tion fraction of the inner stellar medium is likely to be low, which 
would also support the idea that photo-ionisation is the dominant 
source of ionisation. If ionisation is dominated by shocks gener- 
ated by the jet, one would expect that B0108+388 would be rela- 
tively bright in [O III] as its interstellar medium density is appar- 
ently high. B 1358+624 deviates less from the expected relation in 
the right hand panel of Fig.|5|due to its relatively flat X-ray spec- 
trum, which makes that the number flux of UV photons is relatively 
small with respect to the X-ray luminosity. However, since we do 
not know the broad band spectral shape, we do not want to overem- 
phasise this. 

Let us now consider the possible implications of dust absorp- 
tion. The optical depth of dust particles depends on dust particles 
cross sections (ad = ivr 2 with r the physical size of the parti- 
cles, r ~ 0.1 pm) and the column density of dust particles Nd- To 
obtain an order of magnitude estimate we assume that most dust 
particles consist of silicates, and that all silicon is depleted into 
dust. As Nsi ~ 4 x 10~ 4 ./Vh, and a typical dust particle density is 
p — 3.5 g cm -3 , we have Nd ~ 10 _13 A r H, and the optical depth 
should be around adNd — 3 x 10~ 23 -/Vh- This means that dust 
particles absorb an appreciable amount of UV flux if the hydrogen 
column density is comparable to, or exceeds, 10 23 cm -2 , which is 
only the case for B0108+388. We have therefore extrapolated from 
the observed Lx the ionising photon luminosity Nuv using the 

3 Althou gh we do not have kinematic age estimates of PKS1345+125 and 
Mkn668 iGuainazzi et all2004) . reasonable values for the age in fact show 
these galaxies to be too bright in [OIII] compared to the galaxies in our 
sample. In this case the reason is very likely that a large part of the [O III] 
emission is not related to the activity of the central nucleus, as both galaxies 
show evidence of recent merger activity, and are bright infrared sources. 
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Figure S. Left: T he ratio of the [OMT g-^ an d 2-10 keV X-ray luminosity versus the kinematic age of the radio source. The dashed line indicates the average 
value taken from lSambruna et all fi999). Right: Similar to the left hand figure, but instead of Lx the extrapolated UV photon luminosity is used, derived 
from the X-ray luminosity using the power law slope in Tablel3l The arrow indicates approximately how the ratios changes if the UV emission is obscured by 
dust assuming that the dust column density scales with TVh- This would only significantly affect B0108+388. The line shows the expected scaling for a young 
expanding narrow line region (eq.|5j. B1358+624 is indicated with an open symbol, as its age is not based on a direct kinematic age measurement, but on its 
size (TablefTl. 



observed spectral properties. Plotting now I/[oiii] as a function of 
tNuv we see less scatter, certainly if we allow for dust absorption 
(Fig-EJ. In order to bring B0108+388 on the expected relation we 
need a conversion of JVh to dust optical depth that is 16% of the 
above order magnitude estimate. Note that in a log-log plot the ab- 
sorption enters linearly, since \n(Nuv-abs) ~ \n(Nuv) — VdNd- 
Hence, only B0108+388 is likely to be significantly affected by 
dust absorption. 

The large uncertainties in extrapolating from the observed Lx 
to an ionising UV flux makes that we cannot use Fig.|5|(right panel) 
to prove that Eq.[5]is an accurate description, but it makes it at least 
plausible that for GPS/CSO galaxies the [OIII] emission is rela- 
tively low due to an underdeveloped narrow emission line region. 
This is consistent with the idea that GPS-galaxies have AGN that 
switched on around the same time that the radio jets were formed. 

Further support for the idea that GPS/CSO galaxies are in the 
process of creating an extended narrow line region comes from the 
fact that neutral hydrogen apparently extends close to the compact 
radio jets, given the fact that there is strong anti-correlation between 
the ne utral hydrogen column density and jet-size IPihlstrom et alJ 
2003). Note that GPS/CSO galaxies have sizes of the order of a 
few 100 pc, whereas narrow emission line regions can extend up to 
10 kpc. 

3.3 The radio luminosity versus the X-ray luminosity 

The X-ray emission from AGN is thought to come predominantly 
from the immediate vicinity of the central black hole, i.e. thermal 
emission from the accretion disk reprocessed by the hot plasma in 
its vicinity. It is unlikely that synchrotron radiation from the jet 
makes a dominant contribution to the X-ray band. The reason is 
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Figure 6. The ratio of the radio and X-ray luminosity versus the kinematic 
age of the radio source. The dashed line shows th e average ratio for radi o 
loud AGN galaxies with log Lx > 42, taken from Sambru na et al]tl999l) . 
Only sources comparable in X-ray luminosity to our sample were chosen, 
because weaker X-ray sources appear to be relatively radio bright. Like in 
Fig.|5l B 1358+624 is indicated with an open symbol. 
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that, given the typical magnetic fields infered from radio luminosi- 
ties (> 1 mG, IO'Dejll998l) . the synchrotron cooling time rele- 
vant for X-ray synchrotron radiation is in the order of only 2 yr for 
an electron energy of 10 erg. This means that a very small frac- 
tion of the total radio jet volume could produce X-ray synchrotron 
emission. Another potential source of X-ray emission from out- 
side the central region could be inverse Com pton emi ssion by the 
relativistic electron population in the jets (c.f. Belsole et al. 2005). 
Although we cannot totally exclude a significant inverse Compton 
contribution, it seems unlikely to be the case for our sample. The 
reason is that the X-ray emission should in that case come from 
the same locations as the radio emission (the bright regions in the 
jets). However, we would then expect that the measured X-ray ab- 
sorption columns would be more consistent with radio absorption 
measurements toward the jets, which is not the case (section l3~TY 

It is therefore reasonable to assume that, compared to the ra- 
dio and [OIII] luminosity, the X-ray luminosity is more directly 
related to the accretion power of the AGN. Nevertheless, the radio 
and optical luminosities are still indirectly related to the accretion 
power, given the correlatio ns between radio, optical and X-ray lu- 
minosities (e.g. lSambruna et alll99St) . 

It is, therefore, interesting that GPS/CSO galaxies seem on av- 
erage radio bright compared to the X-ray luminosity (Fig.|2j. Given 
the absorption column distribution and low [O III] emission, we can 
ignore the interpretation that the radio emission is relatively bright 
because the interstellar medium is dense in GPS/CSO galaxies, as 
argued by advocates of the "frustrated radio source" scenario. It 
is, therefore, very probable that GPS/CSO galaxies are radio bright 
because they are young. However, within our sample no relation 
between age and radio over X-ray ratio can be seen (Fig.[6}, nor is 
there a correlation with column density. 

Nevertheless, the fact that for the sample as a whole the ra- 
dio to X-ray luminosity is brighter than for other radio-loud AGN 
is at least qualitatively in agreement w ith several radio evolu - 
tion models, such a s Fanti et alj_^995J); iReadhead et al. 1 1996); 
iKaiser & AlexandeJ Jl997l) : lAlexandeJ J2OO0I) and ISnellen et alJ 
(2000). These models describe the evolution of radio jets, with ra- 
dio brightnesses depending on the radial density distribution of the 
interstellar medium. The radio jets are relatively bright as long as 
the they plow through the dense regions of the galaxies, but de- 
cline as soon as they propagate outside the core of the galaxy. The 
difference between the various models is that some assume a den- 
sity distribution described by a King profile ISnellen et al]|200Ct 
Alexander 2000), whereas others assume a density profile falling 
of as power law of radius (Kaiser & Alexander 1997). As a result, 
the non-power law models predict that galaxies in the GPS-phase 
are still increasing in radio luminosity in time, until the jet has 
reached the core radius of ~ 1 kpc, after which the luminosity 
declines. In this case the relatively brightest phase of radio galax- 
ies would be represented by the so-called compact steep spectrum 
sources (CSS), which have more extended radio jets than GPS/CSO 
sources. The power law density models predict that right after the 
jet emerges from the core region the radio emission starts to de- 
cline. 

Our results are inconclusive regarding the details of the early 
evolution of the radio luminosity. However, future X-ray observa- 
tions of CSS galaxies could help to clarify the brightness evolution 
further, as models with a King profile predict that CSS galaxies 
should, on average, have a higher radio to X-ray luminosity ratio 
than GPS/CSO galaxies, whereas models that assume power law 
density profiles predict that CSS galaxies should have a smaller ra- 
dio to X-ray luminosity. 



4 CONCLUSIONS 

We have presented XMM-Newton observations of a sample 
of all the GPS/CSO gal axies with |6| > 20° from the 
Pearson & Readhead ( 1988) catalog, four of which have measured 
kinematic time scales for the jet expansion. All five of the sources 
are detected by XMM-Newton thereby increasing the number of X- 
ray detected GPS/CSO galaxies from 2 to 7. These detections allow 
us to compare the X-ray properties of GPS/CSO galaxies with those 
of other radio loud AGN. 

The results presented here support the hypothesis that 
GPS/CSO galaxies represent the young phases in the evolution of 
radio-loud AGN. The alternative explanation that the radio sources 
are compact due to confinement by an exceptionally high density of 
the interstellar medium in those galaxies, seems extremely unlikely 
in view of the low intrinsic X-ray absorption column densities, 
which ranges from jV H — 4 x 10 21 cm -2 to 6 x 10 23 cm -2 , and 
has a distribution similar to other r adio loud AGN. After submis- 
sion of our manuscript, a preprint by Guainazzi et al. ( 2005) arrived 
at apparently different conclusions based on a sample of five differ- 
ent GPS galaxies observed by Chandra and XMM-Newton. Only 
one of the five GPS galaxies in their sample has Nh< 10 22 cm~ 2 , 
whereas this is 75%±26% for their control sample. Note, however, 
that in our sample three out of five have Nn< 10 22 cm -2 . Taking 
both samples together, this means that four out of ten GPS galaxies, 
or 40%±20% have iV H < 10 22 cm -2 consistent with the control 
sample. 

The fact that the absorption columns toward GPS galaxies are 
consistent with those of other radio galaxies suggests that the in- 
terstellar medium densities in the cores of GPS/CSO galaxies are 
similar to those of oth er radio lou d AGN. A similar conclusion was 
reached by Pihl strom et alJ 12003) based on HI radio absorption ob- 
servations, but the X-ray data provide stronger constraints, as the 
total column density contributes to the absorption, including ion- 
ized regions of the interstellar medium. This may contribute to the 
fact that in all cases the X-ray column densities are higher than the 
HI column densities. A difference between the radio column densi- 
ties and the X-ray column densities is also that the X-ray column is 
measured toward the central source, whereas the HI column density 
is toward the jets, which extend outside the central region. If the dif- 
ference between radio and X-ray column density is dominated by 
those geometrical effects, this would be additional evidence against 
the confinement scenario, since the jets have apparently been able 
to pierce through the dense local regions that contribute most to the 
X-ray absorption column. 

Although the Nh distribution of our sample cannot be dis- 
tinguished from other radio-loud AGN, the ratio of radio to X- 
ray luminosity shows that GPS galaxies have a strong tendency 
to be relatively radio bright. This supports the view that for the 
same thrust of the jets younger radio sources are relatively bright 
iKaiser & AlexandeJ 1997llSnellen et all200d) . However, the data 
is inconclusive concerning whether GPS-galaxies represent the 
most radio-luminous phases in the lifes of radio-loud AGN, as 
would be the case if the source devel ops in a de nsity profi le that 
drops of as power la w with distance iKaiser & AlexandeJ 1 1997t 
|Pihlstr6met'ai] r2003). or whether they are still in their brighten- 
ing phase. This would be the case if the interstellar medium is 
best described by a King profile with a relative uniform density 
with in ~1 kpc of the ce nter and then dropping of as a power law 
(e.g. lSnellenetalj200rj) . In the latter case the brightest evolution- 
ary phase of radio-loud AGN would be represented by the com- 
pact steep spectrum sources (CSS), which have more extended ra- 
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dio emission than GPS/CSO galaxies. A similar study to this one 
concerning CSS-galaxies can clarify this issue. 

Finally, we find that GPS/CSO galaxies are relatively weak 
in [OIII] line emission. This is again in support of the idea that 
GPS/CSO galaxies represent the very earliest stages of the evo- 
lution of radio-loud AGN, since narrow line regions need time to 
build up to their equilibrium size, and young narrow line regions 
are therefore not as bright as fully developed ones. The narrow line 
region is powered by the UV flux of the central source, but also 
shocks induced by the expanding jets are likely to contribute to 
their formation. For those very young radio-loud AGN we advo- 
cate here that their emission line regions are powered by the UV 
flux from the central sources. A case in point is that the relatively 
weakest [O III] source, B0108+388, has also the highest X-ray col- 
umn density, which suggests that a dusty torus is partially blocking 
the UV light. The low [OIII] luminosity therefore shows that the 
birth of the narrow emission line region must coincide more or less 
with the birth of the radio jet. However, we caution that we only 
have a limited knowledge of the nature of the ionization mecha- 
nism for the [O III] line emission, i.e. both shocks from the jets, as 
the UV radiation from the AGN may contribute to the ionization. 
Moreover, compact steep spectrum (CSS) sources, probably repre- 
senting a more advanced evolutionary state of radio galaxies than 
GPS galaxies, show that the forbidden line emission tends to be 
aligned with the radio jet Ide Vries et alTl 999V This phenomenon 
is not well understood, but it should be accounted for if one wants 
to build a more detailed model of the evolution of emission line 
nebulae in radio galaxies. 

In summary, the findings from this X-ray study lends further 
support to the theory that GPS/CSO galaxies represent the early 
phases of radio-loud AGN, in which also the narrow line emission 
nebula is still in the early phases of its evolution. Future X-ray stud- 
ies may help to further clarify the relation between age or jet-size, 
the extent or brightness of the emission line region and the power 
of the X-ray emission. It is important to include also compact steep 
spectrum (CSS) sources in such a study, as they are likely to repre- 
sent the next phase in the evolution of radio-loud AGN. Comparing 
their X-ray to radio luminosity may help clarify whether the radio 
emission declines already during the GPS-phase, or first increases, 
then peaks around the CSS-phase and from then on weakens. 
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